Integrative annotation scores of variants for impact on RNA binding protein activities.

MOTIVATION: The ENCODE project generated a large collection of eCLIP-seq RNA binding protein (RBP) profiling data with accompanying RNA-seq transcriptomes of shRNA knockdown of RBPs. These data could have utility in understanding the functional impact of genetic variants, however their potential has not been fully exploited. We implement INCA (Integrative annotation scores of variants for impact on RBP activities) as a multi-step genetic variant scoring approach that leverages the ENCODE RBP data together with ClinVar and integrates multiple computational approaches to aggregate evidence. RESULTS: INCA evaluates variant impacts on RBP activities by leveraging genotypic differences in cell lines used for eCLIP-seq. We show that INCA provides critical specificity, beyond generic scoring for RBP binding disruption, for candidate variants and their linkage-disequilibrium partners. As a result, it can, on average, augment scoring of 46.2% of the candidate variants beyond generic scoring for RBP binding disruption and aid in variant prioritization for follow-up analysis. AVAILABILITY AND IMPLEMENTATION: INCA is implemented in R and is available at https://github.com/keleslab/INCA.

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model.

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will only be maintained if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74-94% of the variance in aneuploid strains' growth rates is explained by the additive cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of snoRNAs and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

The response to single-gene duplication implicates translation as a key vulnerability in aneuploid yeast.

Aneuploidy produces myriad consequences in health and disease, yet models of the deleterious effects of chromosome amplification are still widely debated. To distinguish the molecular determinants of aneuploidy stress, we measured the effects of duplicating individual genes in cells with varying chromosome duplications, in wild-type cells and cells sensitized to aneuploidy by deletion of RNA-binding protein Ssd1. We identified gene duplications that are nearly neutral in wild-type euploid cells but significantly deleterious in euploids lacking SSD1 or SSD1+ aneuploid cells with different chromosome duplications. Several of the most deleterious genes are linked to translation; in contrast, duplication of other translational regulators, including eI5Fa Hyp2, benefit ssd1Δ aneuploids over controls. Using modeling of aneuploid growth defects, we propose that the deleterious effects of aneuploidy emerge from an interaction between the cumulative burden of many amplified genes on a chromosome and a subset of duplicated genes that become toxic in that context. Our results suggest that the mechanism behind their toxicity is linked to a key vulnerability in translation in aneuploid cells. These findings provide a perspective on the dual impact of individual genes and overall genomic burden, offering new avenues for understanding aneuploidy and its cellular consequences.

Advances in S. cerevisiae Engineering for Xylose Fermentation and Biofuel Production: Balancing Growth, Metabolism, and Defense.

Genetically engineering microorganisms to produce chemicals has changed the industrialized world. The budding yeast Saccharomyces cerevisiae is frequently used in industry due to its genetic tractability and unique metabolic capabilities. S. cerevisiae has been engineered to produce novel compounds from diverse sugars found in lignocellulosic biomass, including pentose sugars, like xylose, not recognized by the organism. Engineering high flux toward novel compounds has proved to be more challenging than anticipated since simply introducing pathway components is often not enough. Several studies show that the rewiring of upstream signaling is required to direct products toward pathways of interest, but doing so can diminish stress tolerance, which is important in industrial conditions. As an example of these challenges, we reviewed S. cerevisiae engineering efforts, enabling anaerobic xylose fermentation as a model system and showcasing the regulatory interplay's controlling growth, metabolism, and stress defense. Enabling xylose fermentation in S. cerevisiae requires the introduction of several key metabolic enzymes but also regulatory rewiring of three signaling pathways at the intersection of the growth and stress defense responses: the RAS/PKA, Snf1, and high osmolarity glycerol (HOG) pathways. The current studies reviewed here suggest the modulation of global signaling pathways should be adopted into biorefinery microbial engineering pipelines to increase efficient product yields.

PKA regulatory subunit Bcy1 couples growth, lipid metabolism, and fermentation during anaerobic xylose growth in Saccharomyces cerevisiae.

Organisms have evolved elaborate physiological pathways that regulate growth, proliferation, metabolism, and stress response. These pathways must be properly coordinated to elicit the appropriate response to an ever-changing environment. While individual pathways have been well studied in a variety of model systems, there remains much to uncover about how pathways are integrated to produce systemic changes in a cell, especially in dynamic conditions. We previously showed that deletion of Protein Kinase A (PKA) regulatory subunit BCY1 can decouple growth and metabolism in Saccharomyces cerevisiae engineered for anaerobic xylose fermentation, allowing for robust fermentation in the absence of division. This provides an opportunity to understand how PKA signaling normally coordinates these processes. Here, we integrated transcriptomic, lipidomic, and phospho-proteomic responses upon a glucose to xylose shift across a series of strains with different genetic mutations promoting either coupled or decoupled xylose-dependent growth and metabolism. Together, results suggested that defects in lipid homeostasis limit growth in the bcy1Δ strain despite robust metabolism. To further understand this mechanism, we performed adaptive laboratory evolutions to re-evolve coupled growth and metabolism in the bcy1Δ parental strain. The evolved strain harbored mutations in PKA subunit TPK1 and lipid regulator OPI1, among other genes, and evolved changes in lipid profiles and gene expression. Deletion of the evolved opi1 gene partially reverted the strain's phenotype to the bcy1Δ parent, with reduced growth and robust xylose fermentation. We suggest several models for how cells coordinate growth, metabolism, and other responses in budding yeast and how restructuring these processes enables anaerobic xylose utilization.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast.

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast.

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

The IV International Symposium on Fungal Stress and the XIII International Fungal Biology Conference.

For the first time, the International Symposium on Fungal Stress was joined by the XIII International Fungal Biology Conference. The International Symposium on Fungal Stress (ISFUS), always held in Brazil, is now in its fourth edition, as an event of recognized quality in the international community of mycological research. The event held in São José dos Campos, SP, Brazil, in September 2022, featured 33 renowned speakers from 12 countries, including: Austria, Brazil, France, Germany, Ghana, Hungary, México, Pakistan, Spain, Slovenia, USA, and UK. In addition to the scientific contribution of the event in bringing together national and international researchers and their work in a strategic area, it helps maintain and strengthen international cooperation for scientific development in Brazil.

Modeling single-cell phenotypes links yeast stress acclimation to transcriptional repression and pre-stress cellular states.

Stress defense and cell growth are inversely related in bulk culture analyses; however, these studies miss substantial cell-to-cell heterogeneity, thus obscuring true phenotypic relationships. Here, we devised a microfluidics system to characterize multiple phenotypes in single yeast cells over time before, during, and after salt stress. The system measured cell and colony size, growth rate, and cell-cycle phase along with nuclear trans-localization of two transcription factors: stress-activated Msn2 that regulates defense genes and Dot6 that represses ribosome biogenesis genes during an active stress response. By tracking cells dynamically, we discovered unexpected discordance between Msn2 and Dot6 behavior that revealed subpopulations of cells with distinct growth properties. Surprisingly, post-stress growth recovery was positively corelated with activation of the Dot6 repressor. In contrast, cells lacking Dot6 displayed slower growth acclimation, even though they grow normally in the absence of stress. We show that wild-type cells with a larger Dot6 response display faster production of Msn2-regulated Ctt1 protein, separable from the contribution of Msn2. These results are consistent with the model that transcriptional repression during acute stress in yeast provides a protective response, likely by redirecting translational capacity to induced transcripts.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions.

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Comparative functional genomics identifies an iron-limited bottleneck in a Saccharomyces cerevisiae strain with a cytosolic-localized isobutanol pathway.

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway re-localization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADH-dependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe-S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

Breadth and Specificity in Pleiotropic Protein Kinase A Activity and Environmental Responses.

Protein Kinase A (PKA) is an essential kinase that is conserved across eukaryotes and plays fundamental roles in a wide range of organismal processes, including growth control, learning and memory, cardiovascular health, and development. PKA mediates these responses through the direct phosphorylation of hundreds of proteins-however, which proteins are phosphorylated can vary widely across cell types and environmental cues, even within the same organism. A major question is how cells enact specificity and precision in PKA activity to mount the proper response, especially during environmental changes in which only a subset of PKA-controlled processes must respond. Research over the years has uncovered multiple strategies that cells use to modulate PKA activity and specificity. This review highlights recent advances in our understanding of PKA signaling control including subcellular targeting, phase separation, feedback control, and standing waves of allosteric regulation. We discuss how the complex inputs and outputs to the PKA network simultaneously pose challenges and solutions in signaling integration and insulation. PKA serves as a model for how the same regulatory factors can serve broad pleiotropic functions but maintain specificity in localized control.

Crabtree/Warburg-like aerobic xylose fermentation by engineered Saccharomyces cerevisiae.

Bottlenecks in the efficient conversion of xylose into cost-effective biofuels have limited the widespread use of plant lignocellulose as a renewable feedstock. The yeast Saccharomyces cerevisiae ferments glucose into ethanol with such high metabolic flux that it ferments high concentrations of glucose aerobically, a trait called the Crabtree/Warburg Effect. In contrast to glucose, most engineered S. cerevisiae strains do not ferment xylose at economically viable rates and yields, and they require respiration to achieve sufficient xylose metabolic flux and energy return for growth aerobically. Here, we evolved respiration-deficient S. cerevisiae strains that can grow on and ferment xylose to ethanol aerobically, a trait analogous to the Crabtree/Warburg Effect for glucose. Through genome sequence comparisons and directed engineering, we determined that duplications of genes encoding engineered xylose metabolism enzymes, as well as TKL1, a gene encoding a transketolase in the pentose phosphate pathway, were the causative genetic changes for the evolved phenotype. Reengineered duplications of these enzymes, in combination with deletion mutations in HOG1, ISU1, GRE3, and IRA2, increased the rates of aerobic and anaerobic xylose fermentation. Importantly, we found that these genetic modifications function in another genetic background and increase the rate and yield of xylose-to-ethanol conversion in industrially relevant switchgrass hydrolysate, indicating that these specific genetic modifications may enable the sustainable production of industrial biofuels from yeast. We propose a model for how key regulatory mutations prime yeast for aerobic xylose fermentation by lowering the threshold for overflow metabolism, allowing mutations to increase xylose flux and to redirect it into fermentation products.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories.

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

The role of stress-activated RNA-protein granules in surviving adversity.

Severe environmental stress can trigger a plethora of physiological changes and, in the process, significant cytoplasmic reorganization. Stress-activated RNA-protein granules have been implicated in this cellular overhaul by sequestering pre-existing mRNAs and influencing their fates during and after stress acclimation. While the composition and dynamics of stress-activated granule formation has been well studied, their function and impact on RNA-cargo has remained murky. Several recent studies challenge the view that these granules degrade and silence mRNAs present at the onset of stress and instead suggest new roles for these structures in mRNA storage, transit, and inheritance. Here we discuss recent evidence for revised models of stress-activated granule functions and the role of these granules in stress survival and recovery.

Genetic variation in aneuploidy prevalence and tolerance across Saccharomyces cerevisiae lineages.

Individuals carrying an aberrant number of chromosomes can vary widely in their expression of aneuploidy phenotypes. A major unanswered question is the degree to which an individual's genetic makeup influences its tolerance of karyotypic imbalance. Here we investigated within-species variation in aneuploidy prevalence and tolerance, using Saccharomyces cerevisiae as a model for eukaryotic biology. We analyzed genotypic and phenotypic variation recently published for over 1,000 S. cerevisiae strains spanning dozens of genetically defined clades and ecological associations. Our results show that the prevalence of chromosome gain and loss varies by clade and can be better explained by differences in genetic background than ecology. The relationships between lineages with high aneuploidy frequencies suggest that increased aneuploidy prevalence emerged multiple times in S. cerevisiae evolution. Separate from aneuploidy prevalence, analyzing growth phenotypes revealed that some genetic backgrounds-such as the European Wine lineage-show fitness costs in aneuploids compared to euploids, whereas other clades with high aneuploidy frequencies show little evidence of major deleterious effects. Our analysis confirms that chromosome gain can produce phenotypic benefits, which could influence evolutionary trajectories. These results have important implications for understanding genetic variation in aneuploidy prevalence in health, disease, and evolution.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations.

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses.

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. These results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

The genetic basis of aneuploidy tolerance in wild yeast.

Aneuploidy is highly detrimental during development yet common in cancers and pathogenic fungi - what gives rise to differences in aneuploidy tolerance remains unclear. We previously showed that wild isolates of Saccharomyces cerevisiae tolerate chromosome amplification while laboratory strains used as a model for aneuploid syndromes do not. Here, we mapped the genetic basis to Ssd1, an RNA-binding translational regulator that is functional in wild aneuploids but defective in laboratory strain W303. Loss of SSD1 recapitulates myriad aneuploidy signatures previously taken as eukaryotic responses. We show that aneuploidy tolerance is enabled via a role for Ssd1 in mitochondrial physiology, including binding and regulating nuclear-encoded mitochondrial mRNAs, coupled with a role in mitigating proteostasis stress. Recapitulating ssd1Δ defects with combinatorial drug treatment selectively blocked proliferation of wild-type aneuploids compared to euploids. Our work adds to elegant studies in the sensitized laboratory strain to present a mechanistic understanding of eukaryotic aneuploidy tolerance.

PKA and HOG signaling contribute separable roles to anaerobic xylose fermentation in yeast engineered for biofuel production.

Lignocellulosic biomass offers a sustainable source for biofuel production that does not compete with food-based cropping systems. Importantly, two critical bottlenecks prevent economic adoption: many industrially relevant microorganisms cannot ferment pentose sugars prevalent in lignocellulosic medium, leaving a significant amount of carbon unutilized. Furthermore, chemical biomass pretreatment required to release fermentable sugars generates a variety of toxins, which inhibit microbial growth and metabolism, specifically limiting pentose utilization in engineered strains. Here we dissected genetic determinants of anaerobic xylose fermentation and stress tolerance in chemically pretreated corn stover biomass, called hydrolysate. We previously revealed that loss-of-function mutations in the stress-responsive MAP kinase HOG1 and negative regulator of the RAS/Protein Kinase A (PKA) pathway, IRA2, enhances anaerobic xylose fermentation. However, these mutations likely reduce cells' ability to tolerate the toxins present in lignocellulosic hydrolysate, making the strain especially vulnerable to it. We tested the contributions of Hog1 and PKA signaling via IRA2 or PKA negative regulatory subunit BCY1 to metabolism, growth, and stress tolerance in corn stover hydrolysate and laboratory medium with mixed sugars. We found mutations causing upregulated PKA activity increase growth rate and glucose consumption in various media but do not have a specific impact on xylose fermentation. In contrast, mutation of HOG1 specifically increased xylose usage. We hypothesized improving stress tolerance would enhance the rate of xylose consumption in hydrolysate. Surprisingly, increasing stress tolerance did not augment xylose fermentation in lignocellulosic medium in this strain background, suggesting other mechanisms besides cellular stress limit this strain's ability for anaerobic xylose fermentation in hydrolysate.